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Abstract 

We study the intermittent fluorescence of a single molecule, jumping from the "light 
on" to the "light off' state, as a Poisson process modulated by a fluctuating envi- 
ronment. We show that the quasi-periodic and quasi-deterministic environmental 
fluctuations make the distribution of the times of sojourn in the "light off' state 
depart from the exponential form, and that their succession in time mirrors environ- 
mental dynamics. As an illustration, we discuss some recent experimental results, 
where the environmental fluctuations depend on enzymatic activity. 
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1 Introduction 

In the last few years there has been an increasing interest for single-molecule 
and, more in general, single-system spectroscopy [1,2J. Usually these systems 
are characterized by intermittent fluorescence, namely, by sudden jumps be- 
tween two states, a light on and light off state. The Probability Density Func- 
tion (PDF) of the times of sojourn or Waiting Times (WTs) r in these two 
states are denoted by the symbols ipon{T) and i'offi,'^)- In the case of Blinking 
Quantum Dots (BQDs) [3ll^ll5]|6] there is a general agreement that these are 
renewal processes, namely, that the jumps between the two states have the 
effect of resetting to zero the system's memory. The PDFs ipon{T) and ipoff{'T) 
are found to be inverse power laws, whose indexes are, according to Ref. [7|, 
Hon = 1-7 and Hoff = 1.75. 

The authors of Ref. |E] afford directions on how to distinguish a Renewal non 
Poisson Process (RNPP), which is a homogeneous in time process from the 
apparently equivalent condition generated by a Non-Homogeneous Poisson 
Process (NHPP), i.e., with a time dependent rate of event productions q{t). 
Let us recall that a (homogeneous) Poisson process is a particular renewal pro- 
cess where the WT-PDF is an exponential function whose decay rate coincides 
with the (constant) rate of event productions. The physical origin of RNPPs 
is still the object of research and debate, whereas that of NHPPs requires 
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however the identification of dynamics responsible for the rate time depen- 
dence. In this Letter we consider the recent work of Ref. [9], which estabhshes 
the influence of the catalytic activity of a fluctuating enzyme on the intermit- 
tent fluorescence of the fluorogenic substrate. In this case the function iponi'T) 
is exponential, and the complexity of the system is signaled by the function 
iljoff{T), which turns out to be a stretched exponential 



with a = 0.15. We argue that this is a NHPP, namely, an example of the second 
class of processes. According to the earlier work [SlllOfllj a non-exponential 
function can be expressed as the sum of infinitely many exponential functions 
p!2llT3] . In Refs. [SflOfll] this condition was realized by means of the stochastic 
modulation of the event production rate q{t). In the present case, we adopt 
a different kind of modulation, which is quasi-periodic, and this allows us to 
establish an evident connection with the popular Stochastic Resonance (SR) 
effect [Hj . SR is a phenomenon produced by making the event production rate 
depend harmonically on time, a condition generating a periodic WT reorder- 
ing: the length of the WTs as a function of the time order reproduces the rate 
time dependence [15]. In this Letter the quasi-periodic modulation has an in- 
ternal origin and is based on the enzyme's conformation changes. The ultimate 
theoretical justification for the adoption of the periodic assumption is given 
by the work of Winfree [T6j, who depicts the fundamental biological processes 
as operating through cyclical processes. Thus, we assume that the enzyme's 
conformation changes are periodic, although non-harmonic, functions of time. 
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2 Slow Modulation 



We assume that the time duration of the on state is negUgible, so that the 
system, after jumping from the off to the on state, jumps back to the ojf 
state almost immediately. Thus, the on state is a kind of flashing snapshot, 
and the distance between two consecutive snapshots is described by i^offi^)- 
Let us omit for simplicity's sake the subscript off and let us assign to ipij) 
the following form: 



V'M^ j n(g)V',(T)dg = V'oe"W , 



(2) 



being V'g(r) — qe~'^'^ the exponential WT-PDF related to the Poisson pro- 
cess with event production rate q and ■0o = 0(/{tor{l/a)) the normalization 
constant. 

The quantity n(g) represents the statistical weight of each ipq{T). We imagine 
q{t) to be a slow function of time, so as to produce many Poissonian events 
with the same rate q. This is called slow modulation condition, and it makes 
it possible to define n(g) on the basis of the enzyme's fluctuations. Each en- 
zyme conformation yields a given q and the slower the enzyme's conformation 
change, the larger the statistical weight of this rate. It is thus obvious that in 
the slow modulation limit the statistical weight n(g) is inversely proportional 
to the modulus of dq/dt: 

dq 



dt 



TIi{q) 



(3) 



The parameter T defines a time scale, driving the slope of q(t) , which can be 
properly adjusted to fit the experimental results. In order to derive a periodic 
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behavior for q{t), we assume the rate q{t) to increase in time from a minimum 
value go to a maximum and very large value g/, from which it decays slowly to 
the minimum value again, and so on, thereby confining the stretched exponen- 
tial distribution to a limited time range, which can, however, include several 
decades, depending on the choice of the parameters qo, qf and T. The switches 
between increasing and decreasing patterns of q{t) correspond to abrupt slope 
changes at times nT, which make q{t) piecewise smooth, and exactly periodic 
as well. However, once a WT is selected we leave it to stretch beyond the ex- 
act time period T, thereby making the exactly periodic slope change produce 
random effects and quasi-periodicity. 

According to Eq. ([2]), iP{t) is the Laplace transform of the function /(g) = 
gn(g). Thus, from the short-time limit, where the stretched exponential reads 

z^(r)~z^o(l-(r/to)"), (4) 

using the Tauberian theorem, we find for n(g), at g 1, the following asymp- 
totic expression 

n(g) = Cg-'^-l (5) 

By plugging Eq. ([5]) within Eq.Q, we obtain the following equations: 
dq g"+^ 



dt TC 



(6) 



The integration of Eq. from t = to t = T yields two different results 
depending on the sign: 
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referring to the cases of increasing (+ sign) and decreasing (— sign) rate, 
respectively. The times Tmc and Tdec are defined by 



T 

^^"'^^i-(go/g/)"+^' 



(g//go)"+^ 



and the parameter r, determining the coefficient a of the stretched exponential 
is given by 

r = ^. (12) 
a + 1 



These analytical results correspond to the correct derivation of the stretched 
exponential under the limiting condition of very slow modulation, namely, 
under the condition that T is large enough. A periodic behavior of q{t) is 
defined by applying suitable time-shifts of the kind t ^ t + nT to the rates 
given in Eqs. ([7]) and (jS]). 



3 Numerical results 



The purpose of this section is to assess by means of artificial sequences how 
large T must be for the WT-PDF ^/'(r) to fit the stretched exponential form 
of Eq. ([1]). On the other hand, we note that the experimental results of Ref. 
[9] do not yield an exactly periodic time ordering. Thus we plan to establish if 
it is possible to find values of T that make ipij) fit the stretched exponential 
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Fig. 1. The WT-PDF V(t) as a function of r. (a) T = 0.1, (b) T = 1, and (c) 
T = 100. The fuh hne is the stretched exponential of Eq. ([T]), with a = 0.2 and to 
determined by a fitting procedure yielding: (a) to = 3.1 ■ 10^^, (b) to = 4.9 • 10^^, 
and (c) to = 5.1 • 10"*^. 

condition without producing a too markedly periodic time ordering, which 
would depart qualitatively from the experimental results of Ref. [9]. 

The numerical algorithm is based on the jump probability in the time interval 
/ = [t,t + dt]: p{t) = q{t)dt, dt being small enough to satisfy the conditions 
p{t) <^ 1 and q{t) nearly constant on the interval /. We used the following 
values of the parameters: a = 0.2, go = 10 and qj = 10^. Thanks to Eq. ([2]) 
the only fitting parameter remains to- Figures [T](a) to[T]^c) show the WT-PDF 
iP{t) for different values of the parameter T. Fig. [TJ^a) refers to T = 0.1. In 
this case the slow modulation condition is violated. This is the case where the 
process can be classified as renewal. 

We see that for both Fig. [H^b) and Fig. [I^c), referring to T = 1 and T = 100, 
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respectively, the agreement between the numerical result and the stretched 
exponential prediction of Eq. ([T]) is accurate throughout the whole time regime. 
Both conditions, therefore, are in principle compatible with the experiment of 
Ref. [9]. To select which of the two conditions is closer to the experiment 
of Ref. [9], let us plot the WTs Tj according to the chronological order, i.e., 
with respect to their order of appearance i. The results are illustrated in 
Fig. [21 for different values of the parameter T. We make the comparison among 
different values of T with two distinct criteria: in the left column we adopt 
different numbers of WTs, namely we refer to different chronological times, 
so as to produce the same number of abrupt slope changes of g(t); in the 
right column, on the contrary, the different plots contain the same number 
of WTs, namely range from 1 to the same maximum chronological time. The 
results of the right column allow us to establish a qualitative comparison with 
the experimental results of Ref. [9]. Furthermore, we see that in accordance 
with the stochastic resonance phenomenon [15], increasing T has the effect 
of making more evident the resonant character of this effect. We see that 
the condition corresponding to T = 1, where the resonant effect is much less 
marked than in the case T = 100, produces the best qualitative agreement 
with the experimental result of Ref. [9]. It is worth remarking that in this 
condition the WT-PDF V^(r) still retains a markedly stretched exponential 
behavior so as to ensure a fully qualitative agreement with the experimental 
results of Ref. |.9.J- 
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Fig. 2. A comparison of the chronological order of WTs for different values of the 
parameter T. On the left side we arranged the same number of peaks (fixed number 
of enzyme fluctuations), while on the right side we compare the number of peaks 
per fixed number of WTs (or system's jumps). T = 0.1, 1, 100 from the top to the 
bottom. 



4 Aging Experiment and Conclusion 



To draw final conclusions we apply to the artificial sequences {U} of Section [3] 
the Aging Experiment (AE) described in the earlier publication of Ref. 
Here we briefiy recall that the method is based on the computation of the 
survival probabilities '^expi'Tita) and \E'ren(T, ^a)- Both survival probabilities 
are aged and derived from time sequence by delaying the beginning of the 
waiting process for an event until a distance ta from each event of the se- 
quence, "^exp is derived by applying this procedure to the time sequence, and 
'^ren is derived from the experimental histogram by assuming that the renewal 
condition holds true. Upon increase of ta, both expiT^ta) and ^E'renl'T, ^a) can 
become much slower than ^q{t) = ^(r, = 0). If ^E'exp(T, ^a) = ^ren(T, ta), 
and both are slower than \E'o(t), this is the sign that the process is renewal 
and non-Poissonian, insofar as in the Poisson renewal process does not 

yield aging. In Fig. [3] and Fig. HI we show a pair of dotted curves, illustrat- 
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ing "^expi^jta) at two different ages, and a pair of dashed curves, illustrating 
^ren('i", ^a) for the Corresponding ages. In both cases the two curves of the 
latter pair generate a decay significantly slower than \1/o(t). This is only a 
consequence of the fact that \E'o(t) is a stretched exponential. We see, how- 
ever, that in Fig. [3] both dotted curves are slower than \1'o(t), although faster 
than the corresponding dashed curves, thereby indicating that a form of ag- 
ing may exist, albeit significantly reduced. This is an indication that many 
renewal events are present [10] . In the case of Fig. IH on the contrary, the two 
dotted curves virtually coincide with "^oi^r), if we neglect the short-time region, 
thereby implying that the deviation from exponential relaxation is caused by 
a so slow modulation as to cancel all the renewal events that may be non- 
Poissonian. We note that the condition of reduced aging of Fig. [3] corresponds 
to T = 1, which, as shown by Fig. [21 generates for the chronological order a 
structure qualitatively similar to that of the experimental results of Ref. [9]. 

We thus conclude that the sequence of WTs emerging from the experimen- 
tal results of Ref. [9] may be a NHPP, generated by a stochastic resonance 
produced by the enzyme's conformation changes. The debate on whether non- 
exponential relaxation is due to inhomogeneities or cooperative molecular mo- 
tion is of wide interest in biophysics (see, for instance Ref. [H]). We note that 
the recent Letter of Ref. [TS] confirms the stretched nature of the WT-PDF of 
Ref. [9j while demanding further theoretical progress to understand the corre- 
lation among events. Then we notice that experimental evidence for a periodic 
behavior of intermittent fluorescence can also be found in [191 ■ think that 
the theory developed in this Letter, along the lines of the earlier work of [15], 
may turn out to be useful to relate the experimental results to the enzyme 
dynamics. 
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Fig. 3. Results of the renewal aging analysis on the survival probability, for aging 
times ta = 0.1,6.1 from the bottom to the top, and parameter T = 1. The aging 
effect is reduced. 




Fig. 4. Results of the renewal aging analysis on the survival probability, for aging 
times ta = 1,10 from the bottom to the top, and parameter T = 100. The renewal 
aging is almost absent in this case. 
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